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Spinless fermions with Coulomb interaction in a square disordered lattice form a new state of the 
matter which is nor a Fermi glass, neither a Wigner crystal for intermediate Coulomb interaction. 
From a numerical study of small clusters, we find that this new state occurs between the two critical 
carrier densities where insulator-metal-insulator transitions of a hole gas have been observed in 
GaAs. 
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PACS: 71.30.+h, 72.15.Rn 

An important parameter for a system of charged par- 
ticles is the Coulomb energy to Fermi energy ratio r s . In 
a disordered two-dimensional system, the ground state is 
obvious in two limits. For large r s , the charges form 
a pinned Wigner crystal, the Coulomb repulsion be- 
ing dominant over the kinetic energy and the disorder. 
For small r s , the interaction becomes negligible and the 
ground state is a Fermi glass with localized one electron 
states, in agreement with the scaling theory of Anderson 
localization. There is no theory for intermediate r s , while 
many transport measurements following the pioneering 
works of Kravchenko et al jjj and made with electron 
and hole gases give evidence of an intermediate metal- 
lic phase in two dimensions, observed Q for instance 
when 6 < r s < 9 for a hole gas in GaAs heteostructures. 
Our study of spinless fermions with Coulomb repulsion 
in small disordered 2d clusters confirms that there is a 
new ground state for those values of r s . In a given clus- 
ter, as we turn on the interaction, the Fermi ground state 
can be followed from r s = up to a first level crossing 
at rf 



A second crossing occurs at a larger threshold 
after which the ground state can be followed to the 



oo. There is then an intermediate state be- 



limit r 

tween rf and rf . In small clusters, the location of the 
crossings depends on the considered random potentials, 
but a study over the statistical ensemble of the currents 
supported by the ground state gives us two well defined 
values for rf and rf: Mapping the system on a torus 
threaded by an Aharonov-Bohm flux, we denote respec- 
tively // and It the total longitudinal (direction enclosing 
the flux) and transverse parts of the driven current. One 
finds for their typical values \It/Ii\ ~ exp — (r s /rf) and 
/; w exp— (r s /rf) with rf < rf . For the Fermi glass, 
the flux gives rise to a glass of local currents and the 
sign of /; can be diamagnetic or paramagnetic, depend- 
ing on the random potentials. For the intermediate state 
(rf < r s < rf ): the transverse current is suppressed 
while a plastic flow of longitudinal currents persists up 
to rf , where charge crystallization occurs. The sign of 
Ii can be paramagnetic or diamagnetic depending on the 
filling factor (as for the Wigner crystal), but does not de- 
pend on the random potentials (in contrast to the Fermi 



glass) . This suggests that a theorem || giving the sign of 
the current in lei could be extended in 2d when r s > rf . 
For a disorder yielding Anderson localization inside the 
clusters, one finds rf ~ 5.8 and rf w 12.4 in agreement 
with the values given by transport measurements which 
we shortly review. 

In exceptionally clean GaAs/AlGaAs heterostructures, 
an insulator-metal transition (IMT) of a hole gas re- 
sults p[] from an increase of the hole density induced by 
a gate. This occurs at r s 35, in close agreement to 
rf w 37, where charge crystalization takes place accord- 
ing to Monte Carlo calculations ||, and makes highly 
plausible that the observed IMT comes from the quan- 
tum melting of a pinned Wigner crystal. The values of 
r s where an IMT has been previously seen in various 
systems (Si-Mosfet, Si-Ge, GaAS) are given in Ref. [Q, 
corresponding to different degrees of disorder (measured 
by the elastic scattering time r). Those r s drop quickly 
from 35 to a constant value r s w 8 — 10 when r becomes 
smaller. This is again compatible with rf « 7.5 given 
by Monte Carlo calculations || for a solid-fluid transi- 
tion in presence of disorder. If the observed IMT are due 
to interactions, it might be expected that this metallic 
phase will cease to exist as the carrier density is fur- 
ther increased. This is indeed the case || for a hole gas 
in GaAs heterostructures at r s w 6 where an insulating 
state appears, characteristic of a Fermi glass with weak 
electron-electron interactions. A similar observation has 
been also reported for electrons in Si-Mosfet at a lower 
r s w 2. The second transition towards Fermi glass can 
only be easily seen if the disorder is sufficient for driv- 
ing localization effects inside the phase breaking length, 
a condition which weakens the metallic behavior of the 
intermediate phase. 

The purpose of this work is to take advantage of exact 
diagonalization techniques possible only for small clus- 
ters at low filling factors. One can question the ability 
of this method to give results valid for large systems. 
In the clean limit, Pikus and Efros B have obtained 
rf w 35 diagonalizing 6x6 clusters with 6 particles, 



close to 



,.w 



37 obtained by Tanatar and Ceperley 



for the thermodynamic limit. This is a first reason to 
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study how one goes from the Fermi glass towards the 
Wigner crystal in small clusters. A second justification 
will be given a posteriori by the agreement between our 
numerical results and the experiments. We consider a 
two-dimensional model of Coulomb interacting spinlcss 
fcrmions in a random potential. It is defined on a square 
lattice with L 2 sites occupied by N electrons. The Hamil- 
tonians reads 

H = -* Yl C * C J + Y Vi7li + T ~~ ( n * n i)> (!) 



<hj> 



2 tr? »* 



where c[ (a) creates (destroys) an electron in the site 
i = (i x ,i y ), J2<ij> denotes a sum restrained to near- 
est neighbors, t is the strength of the hopping terms 
characterizing the particle kinetic energy (t — 1 in the 
following) and is the inter-particle distance for a 2d 
torus (minimum image convention). The random poten- 
tial Vi of the site i with occupation number rii — c\ci is 
taken from a box distribution of width W . The strength 
U of the true physical long-ranged electron-electron in- 
teraction, suitable for insulating phases where screening 
breaks down, gives a Coulomb energy to Fermi energy 
ratio r s = U/(2ty r Ki') for a filling factor v = N/L 2 . 
For U — the Hamiltonian (Q) reduces to the Ander- 
son model of localization, while for t = it reduces to 
the Coulomb glass model, which describes classical point 
charges in a random potential. Exact diagonalization 
techniques for large sparse matrices (Lanczos method) 
are used to study a statistical ensemble of small clusters. 
The boundary conditions are always taken periodic in the 
transverse y-direction, and are chosen periodic, antiperi- 
odic, or such that the system becomes a torus enclosing 
an Aharonov-Bohm flux (f> in the longitunal x-direction. 

Fig. ^ exhibits behaviors characteristic of individual 
small clusters (L — 6, N = 4), where disorder (W = 15) 
yields a strong Anderson localization, as a function of r s 
(u = 1/9 and < U < 50). Looking at the low energy 
part of the spectrum, one can see that, as we gradu- 
ally turn on the interaction, classification of the levels 
remains invariant up to first avoided crossings, visible 
for the ground state at rf , where a Landau theory of the 
Fermi glass is certainly no longer possible. Looking at 
the electronic density pi = (^ol^il^o) of the ground state 
| "Jo) j we have noticed that it is mainly maximum in the 
minima of the site potentials for the Fermi glass, with re- 
strictions coming from kinetic terms and Pauli principle. 
After the second avoided crossing at rY , Pi is negliglible 
except for four sites forming a lattice of charges as close 
as possible to the Wigner crystal triangular network in 
the imposed square lattice. The degeneracy of the crys- 
tal is removed by the disorder, the array being pinned in 
4 sites of favorable energies. The pi of the intermediate 
ground state give the impression of a mixed phase with 
crytalline domains co-existing with glassy domains where 
fcrmions occupy the sites of lowest energies. 
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FIG. 1. Top: Low energy spectrum of a small single clus- 
ter (JV = 4, L = 6, W = 15) (a 1.9r s term has been sub- 
stracted); Bottom left: jumps of 7 at the second crossing 
where the ground state (filled diamonds) and the first ex- 
cited state (empty diamonds) are interchanged. Bottom right: 
Ground state sensitivity D (left scale, empty circles) and num- 
ber of occupied sites £ 3 (rigth scale, filled circles). 



For the same cluster, we have calculated the density- 
density correlation function C(r) = N^ 1 ^2jPiPi-r and 
the parameter 7 used by Pikus and Efros H for char- 
acterizing the melting of the crystal. The parameter 
7 = max r C(r) — min r C(r) calculated for the ground 
state and the first excited state around allows us to 
identify the second crossing with the melting of the crys- 
tal, since 7 = 1 for a crystal and for a liquid. Moreover, 
one can see that the Wigner crystal becomes unstable in 
the intermediate phase, while the ground state is related 
to the first excitation of the Wigner crystal above rY 
(Fig. [I] bottom left). The general picture is reminiscent 
to those found ||] in one dimensional models. As we in- 
crease r s , there are level crossings associated to charge 
reorganizations of the ground state. As in ref. ||, those 
charge reorganizations are accompanied by noticeable de- 
localization effects of the ground state (Fig. [I] bottom 
right). This is shown by a sharp enhancement of the 
ground state sensitivity D = E(0) — E(tt) measuring the 
change of the ground state energy when the boundary 
conditions are twisted in the x-direction, i.e. more fun- 
damentally the ability of the ground state to a support 
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a persistent current when the system forms a torus en- 
closing a flux. The second quantity demonstrating the 
derealization effect at the crossings is the participation 
ratio £ s — N 2 (J2i Pi)^ 1 of the ground state, i.e. the 
number of sites that it occupies. Fig|l] is representative 
of the ensemble, with the restriction that the location of 
the crossings fluctuates from one sample to another, as 
observed || in Id, as well as the sign (paramagnetic or 
diamagnetic) of D below rf , in contrast to Id. 

Keeping the boundary condition periodic in the trans- 
verse y-direction and such that the system encloses a flux 
cj) (in dimensionless units) in the longitudinal x-direction 
(c/> = 7r corcsponds to anti-periodic condition), one drives 
a persistent current of total longitudinal and transverse 
components given by: 



h = - 



dE{4>) 



yr y. 

and I t = ^ 



respectively. The local current flowing at the site i in the 
x-direction is defined by 



It = 2Im«*„|ci. +1) 
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and by a corresponding expression for If . The response 
is paramagnetic if Ii > and diamagnetic if Ii < 0. 
Legett theorem |3j , based on a simple variational possib- 
lity for the ground state wave function ^o, states that 
the sign of the response is given by the parity of N. 
(— 1) N (E(0) — E(n)) is always positive in Id for all dis- 
order and interaction strength. The proof is based on 
the nature of "non symmetry dictated nodal surfaces ", 
which is trivial in f d, but which has a quite complicated 
topology in higher d. Moreover, in 2d and small r s , the 
sign of /; depends on the site potentials. 

We have separately studied the paramagnetic and dia- 
magnetic samples of an ensemble of 10 3 clusters with 
again L = 6, N — 4 and W — f5. We obtain lognor- 
mal distributions for all values of r s for a large disorder 
(W = 15), as illustrated by Fig.|[ The dependence on 
r s of the averages and variances of the logarithms of the 
paramagnetic D + = E(0) — E(ir) > and diamagnetic 
D_ = E(0) — E(n) < responses are given in Fig|| The 
log-averages exponentially decay as D_ cx exp— (r s /rf) 
and D + cx exp-(r s /rf ) with rf « 5.8 and rf « 12.4. 
The decay of the typical values are accompanied by an 
increase of their fluctuations. The r s -dependence of the 
log- variances are shown in the insert, showing that the 
variances and the averages are proportional. This re- 
minds us a similar relation B for the log-conductance in 
a 2d Fermi glass at r s = 0. The variances of log|/t/7;| 
and log D + increase as r s /rf , and as r s /rf above rf re- 
spectively. Both the average and the variance of log \D + \ 
remain nearly constant in the intermediate phase. 
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FIG. 2. Distribution of the ground state sensitivity D- of 
the diamagnetic samples of an ensemble of 10 4 clusters at 
r a = 1.7 fitted by a log-normal. 




FIG. 3. Log-averages of D+ cx exp — (r s /12.4) (paramag- 
netic, empty circle) and \D- \ cx exp— (r s /5.8) (diamagnetic, 
filled square) given by the left scale, and of the ratio \h/Ii\ 
(empty diamond, right scale). Insert: variances of log|D_|, 
of log |it //; | cxr a /5.9 and of log|_D+| oc r 8 /11.9. 

A transition from glassy towards plastic flow as r s in- 
creases has been found []l0| by Berkovits and Avishai. It 
is likely that the existence of diamagnetic samples comes 
from disordered arrays of local loops of current which 
should give rise to a transverse current I t . Studying the 
distribution of It and J; at (f> = 7r/2, we find that there is 
indeed a non zero transverse current It when r s is small. 
Its sign is random over the statistical ensemble, and the 



ratio log \It/h 
plotted in Fig. 
actly as D_. 



is normally distributed. The log-average 
i confirms that \h/Ii\ oc exp— (r s /rf ) ex- 



The fact that /; becomes paramagnetic independently 
on the microscopic disorder when r s > rf does not mean 
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that Coulomb repulsions always yield this response. For 
instance, 4x6 clusters with N = 6 particles become al- 
ways diamagnetic when r s increases. One just concludes 
that there is a rule for r s > rf giving the sign of the 
response in 2d, as the one found by Legett in Id, which 
does not depend on the random potential. 
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FIG. 4. Flux dependence of ground state and first excited 
state energies for r a = (filled symbols, left scale) and r s = 42 
(empty symbols, right scale) for a single cluster. A suitable 
energy translation has been done to put the curves in the 
same figure. Insert: ensemble average of j3 = DoD\/\DqD\\. 



to Finkelstein's approach |nj , we underline that the spins 
are neglected in our model, and we cannot confirm that 
the intermediate ground state is a metal for a disorder 
energy to kinetic energy ratio as large as W/t = 15. We 
have mostly seen noticeable derealization effects near 
the crossing points. We cannot exclude that weaker dis- 
orders or spin effects are necessary for having a metal. 
We also note that the possibility of a superconducting 



state has been discussed 14 1, and not entirely ruled out 
in the clean limit |l5| ] . A study of the temperature depen- 
dence of the ground state sensitivity and of the system 
size dependence of the gap excitations at a fixed r s will 
be useful for clarifying this issue. From these small clus- 
ter studies, we conclude that there is a new state of the 
matter, clearly separated from the Fermi glass and from 
the Wigner crystal, identified by a plastic flow of currents 
and a magnetic response with a sign independent on the 
microscopic disorder. The obtained critical r s factors are 
surprisingly close to those given by Ref. Q). However, 
the critical r s can depend on W , as already known for 
, and might have small finite size corrections (a study 
in 8 x 8 clusters with N = 4 gives rf « 6.2 instead 



of D_ 
of 5.8). 
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It is also interesting to study the currents supported 
by the first excited state. In a Fermi glass, if the ground 
state current Ii(Eq) has a certain sign, the current Ii{E\) 
carried by the first excitation has very often the oppo- 
site sign. We have noticed that the opposite behavior 
becomes more likely when r s increases: the currents car- 
ried by the ground state and the first excited state flow 
in the same direction. An illustration is given in Fig. ^, 
where the flux dependence of the two first levels at r s = 
and r s — 42 of the same cluster is shown. The tran- 
sition from anti-correlated towards correlated flux de- 
pendence can be seen on the average of D Di/\D n Di\ 
which precisely changes its sign at , as shown in 
the insert. Fig.|| also confirms that the approximation 
—dEj/dcj) « Dj = Ej(0) — Ej(w) makes sense. 

In summary, we have shown that a simple model of 
spinless fcrmions with Coulomb repulsion in a random 
potential can account for the two critical densities of 
holes where insulator-metal-insulator transitions occur. 
This gives an important confirmation that the inter- 
mediate phase observed in two dimensions results from 
Coulomb long-range repulsions, and makes unlikely sce- 
narios based on spin-orbit scattering |pH or on some spe- 
cial single electron interface properties In contrast 
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